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Abstract 



The spectral functions and light-cone momentum distributions of protons 
and neutrons in 3 He and 3 H are given in terms of the three-nucleon wave 
function for realistic nucleon-nucleon interactions. To reduce computational 
complexity, separable expansions are employed for the nucleon-nucleon po- 
tentials. The results for the light-cone momentum distributions suggest that 
they are not very sensitive to the details of the two-body interaction, as long 
as it has reasonable short-range repulsion. The unpolarised and polarised 
structure functions are examined for both 3 He and 3 H in order to test the 
usefulness of 3 He as a neutron target. It is found that the measurement of 
the spin structure function of polarised 3 H would provide a very clear test of 
the predicted change in the polarised parton distributions of a bound proton. 



I. INTRODUCTION 

It is well known that a polarised 3 He target can be used as a polarised neutron target. 
The question we would like to address is how good a polarised neutron target it is for 
the determination of the neutron spin structure function, gi, in deep inelastic scattering. 
There are two questions that play a central role in resolving this problem. The first is the 
sensitivity of the light-front momentum distribution to the three-nucleon wave function. For 
this we need to calculate the spectral function for realistic tri-nucleon wave functions. The 
second question is a consequence of the fact that the neutron structure function is small 
in comparison with the proton structure function. This raises the question of the accuracy 
with which one can extract the polarised neutron structure function from 3 He. 
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To examine these questions we need first to calculate the three-nucleon wave function 
for a "realistic" nucleon-nucleon potential. To simplify the problem computationally, we 
consider a separable expansion [|IJ of the Paris potential (which we call PEST) , that gives 
the same three-nucleon observables as the the original Paris potential in a full multi-channel 
Faddeev calculation |3],|J. For comparison we consider two other classes of potentials. The 
first is a rank one unitary pole approximation (UPA) || to the Reid Soft core potential || . 
This has the property that it reproduces the position and residue of the poles in the 1 So 
and 3 Si- 3 D! channels - i.e., it reproduces the original potential's deuteron wave function. 
As a result, it incorporates the short range behaviour of the original interaction. The 
second is a Yamaguchi type potential with a D-state probability of 4% and 7% [[/J. These 
potentials do not include the short range repulsion that is commonly present in nucleon- 
nucleon interactions. 

In Sec. 2, we present the procedure used to determine the three-nucleon wave functions 
for these potentials, as well as the corresponding three nucleon observables. By comparing 
the results for these three classes of potential, we are able to determine the importance of 
short range correlations and the contribution of higher partial waves to the neutron and 
proton spectral functions and therefore to the light-cone momentum distributions. Since 
we will be considering both 3 He and 3 H, we have chosen to work in an isospin basis and 
therefore neglect the contribution of the Coulomb interaction to the 3 He wave function. We 
do, however, estimate the effect of neglecting the Coulomb correction on the momentum 
distribution and therefore the structure functions. 

In order to analyse the deep inelastic structure functions of A = 3 nuclei, we need to 
determine the neutron and proton spectral functions. This is detailed in Sec. 3. Here we 
compare the results for various two-body potentials, finding that the light-cone momentum 
distribution is not sensitive to the details of our three-nucleon wave function. In Sec. 4 we 
turn to the structure functions and examine the ratio of the structure function in the three- 
nucleon system to that in the deuteron (the EMC effect) for the different interactions. We 
also examine the possible implication of neglecting the Coulomb interaction in 3 He. This 
opens the way for us to study the sensitivity of the unpolarised and polarised structure func- 
tions to the quark distributions in the proton and neutron and the possibility of extracting 
the neutron spin structure function from polarised 3 He data. Finally, in Sec. 4 we present 
some concluding remarks. 

II. THE THREE NUCLEON WAVE FUNCTION 

For the three-nucleon problem we can determine the non-relativistic wave function by 
solving the Faddeev equations exactly for any realistic two-body interaction. However, to 
simplify the computational aspects of the problem, with no sacrifice in the quality of the 
wave function, we turn to separable expansions that have been extensively tested J§|J. 
This will result in a three-nucleon wave function that can be used to calculate the spectral 
function and the light-cone momentum distribution. In the present section we detail the 
three-nucleon formalism required to evaluate the wave functions for 3 He and 3 H. 
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A. Notation 

With the extensive literature on the Faddeev equations || and their use in the three- 
nucleon system, we restrict ourselves here to a summary of the notation used in the present 
analysis. The Faddeev decomposition of the three-nucleon wave function is given by 

I* > = M + M + M = {e + (123) + (132)} |^ 3 ) • (1) 

Here "e", "(123)" and "(132)" are members of the permutation group of three objects, with 
e being the unit element (i.e. e \ip a ) = \<p a )) and the other two being cyclic permutations of 
{1, 2, 3}. The second equality results from the requirement that we have identical particles, 
the wave function is then invariant under any cyclic permutation of our particles. Since we 
have a system of identical fermions, the total wave function must be antisymmetric under the 
exchange of any two particles in the system. This requirement leads to following conditions 

(a/3) \<p a ) = - \<pp) , 

(a/3) \cpfi) = - \cp a ) , (2) 
(a/3) \<p 7 ) = - |</? 7 ) • 

In the above equations a, j3 and 7 are indices running from 1 to 3, and always different from 
each other, and (a/3) is again a member of the permutation group of three objects which 
exchange particles a and (3 leaving the third one unchanged. Since we are dealing with 
a three-body problem, there will be only two independent momenta in the centre of mass 
frame. All the particles have spin and isospin | and one must account for their orbital angular 
momentum. We briefly summarise the quantum numbers and momenta used throughout 
this paper: 

• N a is a set of quantum numbers describing a three body channel from the point of 
view of the particle a, which is the spectator; the set is unique for each channel. 

• £ a is the orbital angular momentum between particles /3 and 7. 

• L a is the orbital angular momentum between particle a and the centre of mass of the 
system consisting of particles /3 and 7. 

• 3m Jpi h are ^ ne spins of each particle. 

• T a , T/3, ty are the isospins of each particle. 

• p a is the momentum of particle a in the centre of mass frame. 

• q a is the relative momentum of the pair of particles (3 and 7, defined as q a = {p 1 —p/s)/2. 

• / and J are respectively the total isospin and total angular momentum of the system. 
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B. The partial wave expansion 



We now turn to the partial wave expansion of our wave function. To minimise the number 
of coupled Faddeev equations, having truncated the interaction to a set of partial waves, we 
have used the following coupling scheme: 



Jf3+Jj = S a , 

If) ~\~ T"y %(X1 



(■a Set Jai J a Jet ^ai L a -\- S a 
1>a let I ) 



J. 



which is known as the channel coupling scheme. With this coupling scheme the complete 
set of quantum number N a describing a three body channel is; N a = {i a , s a ,j a , S a , L Q }. A 
subset of these quantum number that describe the two-body channels is; n a = {i a , s a , j a }, 
and therefore N a = {n a , S a ,L a }. We have not included £ a in the set of quantum numbers 
since the tensor force mixes values of £ a . This allows us to define the angular momentum 
and isospin basis as 



eti ( Jfii J-y) Set) Jcti Jet ]S a }J) \[(ip, ^ 7 )^ a , i a ] I) 



(3) 



These basis states satisfy the following orthogonality relation (fl JI 



e a N 



QJI 



Sc J a 0, 



We are now in a position to write the partial wave expansion of the total three-nucleon 
wave function as 



i a N a 



u 13 

u e a N a 



(4) 



where 
wave -! 



lA\^ Na J is defined as the radial part of the wave function corresponding to the partial 



C. Separable potential 

To reduce the dimensionality of the Faddeev integral equations from two to one, and in 
this way simplify the three-body wave function, we have employed a separable expansion of 
the nucleon-nucleon interaction. Our potential for the interaction of particles j3 and 7 in a 
given partial wave is of the form || 



~[/~n a 
V p pi 

a 5 a 



9e a / A t a i> a \9i> a 



(5) 



where g^"j is a "form factor" and X^p, is the strength of the potential in that partial wave. 
By taking £ a 7^ i' a we can accommodate a tensor interaction, as in the case of the 3 Si- 3 D! 
nucleon-nucleon channel. The above expression for the potential is for a rank one potential. 
To incorporate higher rank potentials, we turn the strength X^p, into a matrix and as a 

result ge"^ is a row matrix. In resorting to separable expansions, we have taken the view 
that the expansion is a numerical procedure analogous to the use of quadratures. However, 
a low order expansion, such as the UPA or the use of a separable potential, is justified on the 
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grounds that it generates the same analytic structure in the amplitude (i.e., bound or anti- 
bound state poles) as a corresponding realistic potential. || The use of a separable potential 
gives rise to a separable t-matrix that satisfies the Lippmann-Schwinger (LS) equation; 



t a (E) =V a + V a G (E) t a (E) = (1 - G (E) V a )~ l V a 



(6) 



with Gq(E) = (E — Hq)^ 1 the two-body Green's function. It is simple to show that the 
separable t-matrix in a given partial wave, resulting from a solution of the LS equation, is 
of the form 



(7) 



where the form factor 



g™* ) is identical to that used in the separable potential. The function 



t%™(, (E), in a given channel, can be a written in matrix form as 



(8) 



This separability of the t-matrix will allow us to reduce the dimensionality of the Faddeev 
integral equations from two to one after the partial wave expansion described in Eq. (01). 



D. The three- nucleon wave function 

Having determined the structure of the two-body amplitude, we now turn to the wave 
function for the three-nucleon system. The Schrodinger equation for this system is 

3 

(E - H ) |tf > =V\V) = J2 V <* I*) • ( 9 ) 

a=l 

This can be rewritten in a form that suggests the Faddeev decomposition stated in Eq. ([[]), 
i.e., 

3 3 

|*) = G (E) V |tt> = <T G (E) V a |tt> = £ \cp a ) . (10) 

a=l a=l 

Here, Gq{E) = (E—H G )~ l is the three-body Green's function. We now can write an equation 
for the Faddeev components of the wave function as 

\<p a ) = G {E) V a |*> = G (E) V a \<p a ) + J] G Q (E) V a |^ 7 ) . (11) 

With the help of Eq. (H) , the set of coupled integral equations for the Faddeev components 
of the wave function, \ip a ), becomes 

\<p a ) = G (E)T a (E)(\^) + \^)) . (12) 

Here T a (E) is the t-matrix for particles j3 and 7 in the three-particle Hilbert space, which 
is related to the two-body amplitude considered in the last section by 
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T a (E) =t a (E-e a ) 



(13) 



where e a is the energy of the spectator particle a in the three-body centre of mass.[] 

In Eq. (|i2|) we have a set of coupled integral equations, known as the Faddeev equations, 
for the three-body bound state. For the three-nucleon system, where we have identical 
Fermions, we take advantage of the anti-symmetry, as given in Eq. (H), and the fact that 
(^j)T a = T a (Pj) = —T a , to reduce the Faddeev equations to 



\<p a ) = G Q (E) T a (E) (1 - (/3 7 )) \(p p ) = 2G (E)T a (E) \<p p ) 



(14) 



with a (3. To recast this equation into a form that will admit numerical solutions, we 
need to first partial wave decompose the Faddeev equations and take into consideration the 
separability of the two-body amplitudes. This can all be achieved by partial wave expanding 
the two-body amplitude in three-body Hilbert space in terms of the angular momentum 
states defined in Eq. (|3]) [[[(J 



r a (E)= £ fdp aP l \n^ Na]Pa )t^ a (E-e a ) ( Pa ,n J e / No 

i> e< „ 



'o'a o 



J2 J dpaP, 



t a l' a o 

N a JI 



where e a = ^-p\ and 



qJI . n a 



n JI 



Pa 



We now can write Eq. ( |14|) as 

oo 
l a f a o 

N a JI 

OO 

= 2G (E) £ JdPaA \^V a N a ;9 n c) rt\(E - e a ) Xj^(p e 



VP 



*■ <* 

N a JI 



with the spectator function, X^ e (p a ), satisfying the equation 

Xil ala {p a ) = (g^KN a \pp) 

oo 

= 2E/ dpppl Z^ Na . £ ^(p a , Pl3 ;E) 

o -el „ 



o 



(15) 



(16) 



(17) 



(18) 



1 For the three-nucleon system in a non-relativistic formulation, e a = zlrPa, where m is the nucleon 



mass. 
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where 



z elN a] i l3 N l3 (Pa^Pi3', E) = ^g^;^ Na G (E) SlipNp'iQip) ■> 



.11 



(19) 



with a 7^ p. In Appendix [A] we give an explicit expression for Z^ Na . e N . for the coupling 
scheme used in the present analysis PJTD| . In Eq. ( |TS| ) we have a set of coupled, homogeneous, 
integral equations for the spectator wave function, X^/ aea (p a ), which we can use to construct 
the total wave function. Here, we note that the spectator wave function is only a function of 
the momentum of the spectator particle and the energy of the system, which is the binding 
energy of 3 He or 3 H. We now turn to the total wave function for the three-nucleon system. 
Making use of the orthogonality of the angular functions, Q 
radial wave function, defined in Eq. (ffl), as 



ji 



we can write the total 



where 



U N a i a 



n JI 
n JI 



Pa) + (Q J e*N a VP + V? 7 



n JI2 



(20) 



'h a N 



2G (q a , Pa ; E) g n £(q a ) £ ^t\( E ~ (P«) 



(21) 



-i 



with G (q a ,Pa',E) = E — ^ ( q 2 a + |p^,j . The second component of the radial wave 
function in Eq. (|20|) is given by 



Ve„N„(P<*, Q a ) = \Paq 

\Paqon ^ 



JI 

£ a N a 



p J2 J d Z T iiN a ^Np (p«, p'^ x ) niphp (p|9> ej) 



+i 



h*fi -i 



(22) 



where P = \ 



1 - 



and 



P'p =ql + \pI + q a Pa^ , $ = ~9« + ^pl ~ ^q a Pa£ 



x = ~7p{l Pa + qa ^ ■ (23) 



The function r/J^.^ N is given in Appendix |A]. We only observe here that the expression 
for F{* Na . i/3N differs from that for Z^ Na . t N by the absence of the separable potential form 
factors and the three-body Green's function. The normalisation of the total wave function 
is then given by 
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= 3(<p a \<Pa) +6{<p a \<pp) 



n JI2 

'k a N 



(24) 



Here the sum is restricted by the two-body partial waves included in the Faddeev equations. 
Since the partial wave expansion of the total wave function involves an infinite sum, we need 
to truncate this sum such that the normalisation evaluated by the truncated sum, that is: 



U 



ji 

N a e a 1 1 



(25) 



agrees with the result of Eq. In this way we ensure that our total wave function includes 
all the partial waves dictated by the two-body interaction. 



E. Numerical results 

As a first step in the determination of our wave function, we calculate the binding energy 
of the three-nucleon system for the class of potentials being considered. For the UPA to 
the Reid Soft core and the Yamaguchi potentials the interaction is restricted to the 1 So and 
3 Si- 3 D! channels. This reduces the homogeneous Faddeev equations to five coupled integral 
equations for the spectator wave function. For the PEST potentials the number of coupled 
channels depends on the rank of the interaction in a given channel and the number of partial 
waves included. To get the optimal representation of the Paris potential we need to have 
achieved convergence in the rank. This varies from channel to channel. In all cases the 
rank has been chosen in such a way that the binding energy for a given number of channels 
has converged and is in agreement with the results of calculations using the Paris potential 
directly . In Table | we present the result for the binding energy for the three classes of 
potentials. For the PEST potentials we have taken the 5, 10, and 18 channel potentials. The 
18 channel calculation corresponds to including all nucleon-nucleon channels with J < 2. 
This will allow us to examine the contribution to the spectral function from higher partial 
waves. Here we observe that the Yamaguchi potentials overbind the three-nucleon system, 
while the UPA and PEST potentials underbind. Since the binding energy determines the 
long range part of the wave function, this difference allows us to examine the sensitivity of 
the structure functions to the binding energy and therefore to the tail of the wave function. 
A comparison of the PEST five channel and the UPA suggests that the difference between 
these two models is minimal. In fact, that is the case for most realistic potentials that do 
not include energy dependence. The higher partial waves in the PEST potential seem to 
have a small but significant contribution to the binding energy. Here again, this potential, 
in common with all realistic potentials, tends to under-bind the three nucleon system. The 
solution to this problem may involve the short-range, velocity dependence of the two-nucleon 
force [[□]], as well as a genuine three-body force |T^| . 



Since we have neglected the Coulomb contribution to the energy of 3 He, and our more re- 
alistic potentials under-bind the three nucleon system, we have chosen to adjust the strength 
of the x So interaction to reproduce the experimental binding energy of both 3 He and 3 H. 
This procedure does not effect the deuteron wave function, but could have some influence 
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on the continuum wave function in the 1 So- In this way, we may estimate the error in ne- 
glecting the Coulomb energy for 3 He, and the possible error in the tail of the wave function 
due to underbinding of the three nucleon system. The contribution of this correction will be 
discussed when considering the spectral functions and light-cone momentum distributions. 



III. LIGHT CONE MOMENTUM DISTRIBUTION 

Before we proceed with the discussion of light-cone momentum distributions, we should 
establish the relation between the cross section in charged lepton scattering and the light- 
cone momentum distribution. The cross section for the scattering of a charged lepton with 
a nucleus is proportional to the product of the leptonic tensor L^ u with the hadronic tensor 
W^ v . For an unpolarised hadronic system of spin 1/2 (i.e. free nucleon, 3 He and 3 H) the 



hadronic tensor has the following form [|l3| --|l5|j 



W ^ = \llj d ^ qX (PS\Ux)U0)\PS) (26) 

/(_1, + ^)h, + (p„-^)(p„-^)§, 

where P is the four momentum of the hadronic system, S is its polarisation and M is its 
mass. Here, J is the electromagnetic current, and q the four momentum of the virtual 
photon. Finally, W\ and W 2 are the form factors of the hadronic system. In deep inelastic 
scattering, one prefers to use the structure functions F\ and F 2 instead. The relation between 
the form factors and the structure functions is the following 

Fi = MWt, F 2 = ^~W 2 . (27) 



The leptonic tensor for unpolarised scattering has the following structure [13-fOJ 

L ^ = \Y^ s')ln u ik, s)u(k', s)^ v u(k, s), 

s,s 

= 2 + K k » - 9»uk ■ k') , (28) 
with k (k') and s (s ) the initial (final) four momentum and polarisation of the lepton. 

For polarised scattering one does not average over the initial polarisation and the re- 
sulting tensors then have two parts; a symmetric part, identical to those of Eq. ( |2T)| ) and 
Eq. (|28|) , and a new antisymmetric piece that is related to the polarisation. The antisym- 
metric part of the hadronic tensor contains two new form factors, G\ and G 2 , which are in 
turn linked to two new structure functions, g\ and g 2 . 

The convolution formalism gives a prescription, valid under certain conditions, to link 
structure functions of complex hadronic systems to structure functions of free nucleons 
16| , [T7jl . In this formalism, the nucleon light cone momentum distribution in a nucleus plays 
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a central role, in that it relates the in-medium structure function to the nucleon structure 
function. This relation takes the form of a convolution integral given by (see Ref. ||18|| ) 



M A / \ 

F$(x,C?)=f m dyf(y)F N (*,Q 2 ) . (29) 



Here, F N (F$) is the free (in medium) structure function, / is the nucleon light cone 
momentum distribution, Ma and m are the masses of the nucleus and of the free nucleon 
respectively, finally, x is the traditional Bjorken variable and Q 2 is the momentum transfer 
squared (Q 2 = —q 2 ). The above relation is valid for the leading twist of the structure 
functions, which is why f(y) has no Q 2 dependence. Another important assumption made 
in this formula is the impulse approximation, namely the assumption that the structure 
function of an off-shell nucleon is equal to the structure function of an on-shell nucleon. A 
more complete discussion about problems raised by this assumption can be found in Ref. 

The nucleon light cone momentum distribution in a nucleus, f(y), is the probability to 
find the nucleon in the nucleus with a given fraction of the total momentum y(= p + /P + ) 
of the nucleus on the light front. As a result, one readily see that Eq. (|29|) has a simple 
interpretation. The structure function in the medium is the sum of all possible values of the 
free nucleon structure function, weighted by the probability of finding the nucleon with a 
given momentum fraction y. In this section, we will show how to determine the light cone 
momentum distributions for the neutron or proton in the three-nucleon system. 

Since the light cone momentum distribution is essentially the probability of finding a 
given nucleon with a particular fraction of the momentum of a nucleus, it should be related 
to the spectral function of the nucleon in that nucleus. In the instantaneous frame the 
spectral function is the combined probability of finding a nucleon with a given momentum k 
while the remaining nucleus is in a state A. We denote this spectral function by S\(k). The 
light cone momentum distribution is then a sum over all possible states A, and all possible 
k that are compatible with the fraction of momentum y. This is given by 

M = Y.J d4k ( 1 + ]p) s (y- ^Pj Sx(k) . (so) 



In some cases (see Ref. ||13|| ) a light cone momentum distribution is defined for each state 
A. In Eq. (§§) the factor (1 + k 3 /k°) is called the flux factor. It is a relativistic correction 
arising from the fact that we are using a light front formalism Light cone momentum 

distributions, as well as spectral functions, can also be defined for polarised nucleons. In the 
following section, we will concentrate on the unpolarised spectral function and merely state 
the results for the polarised nucleon spectral function. 

We note that the calculation of the nucleon momentum distributions presented here is 
very similar in spirit to the pioneering work of Ciofi degli Atti and Liuti [ 2TJ . That work used 



a wave function based on variational method, rather than the Faddeev equations. While the 
variational approach is designed to produce an accurate estimate of the binding energy of the 
system, one must work harder to obtain an equally accurate wave function. Indeed, for the 
tri-nucleon system this has led to the necessity to explicitly correct the proton momentum 
distribution, as described in Ref. We are not aware of a similar correction being applied 
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to the neutron momentum distribution. In any case, it appears to us that it is worthwhile to 
make the calculation with a different technique. In addition, we can study the dependence 
on the assumed two-nucleon force explicitly. 



A. The spectral function 



To determine the light cone momentum distribution we need to know how to compute 
the spectral function. For the unpolarised case, the "diagonal spectral function" is given by 



S x {k) 



2J A + 



2J a + 



^— (y,VA\ai tN {k)a a , N {k)\y,(T A ) 6(k°-(m + ex-Tr) 

1 I ^2 

— — - ^ \(<p,<Tb a^ N (k) *,o-x) 5 (k° - (m + e A - Tr A )) . 



(31) 



Here, a a) is the wave function of the initial nucleus A with spin, J4, and spin projection, 
a a, along the z-axis, while |0, 07,) is the wave function of the A — 1 system in the state ov 
The sum over o\, is restricted to those states allowed by the energy conserving 5-function. 
The operator a\ N (k) is the creation operator for a nucleon iV (proton or neutron) with spin 

projection a and momentum k. 

In the following we will note the product N (k) a a ^(k) as the familiar number density 

operator p a ,N{k) and we will define it in a way similar to Ref. For example, the density 
of protons with spin +1/2 along the z-axis and momentum p, (Pp{p)) , in a tri-nucleon, is 
defined by 



with 



(32) 



1 ^ r 

2 E E J ^ (p, q) \pti\ * , ^ (p, g) 



f"A 2=1 



Pp,i 



l + r 3li ) (1 + a 2 



(33) 



In Eq. ( PB| ) one can recognise the number density, in the sense of Ref. The other density 
operators which we may use are 



Pp,i 
Pn,i 
Pn.i 



;i + r 3 ,i)(l-tr z , i ) 

2 2 

;i -r 3|i ) (l + <r z ,j) 

2 2 

^1 - r 3 ,i) (1 - a Zji ) 
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Using the notation of section 2, and more specifically Eq. (£|), we can rewrite Eq. (|32|) in a 
slightly different way, showing explicitly how we conduct this computation with our wave 
function 



(pt(p)) = l E 



2 



E / d2 V { n e! Na , °a (p, §) |p+; | ^//^ > (p, q) ) 
x (| dqq 2 (lAl J aNa (p, g) (ZY//^ (p, g)) 



(34) 



B. The case of 3 He 

3 He is one of simplest nuclei, along with 3 H and deuterium. It consists of 2 protons and 1 
neutron. If we measure the light-cone momentum distribution of the neutron, the remaining 
two protons can only be in a scattering state, since there is no bound state of two protons. 
On the other hand, if we measure the light cone momentum distribution of the proton, the 
remaining two nucleons are a proton and a neutron, which can be in either a bound state, the 
deuteron, or a scattering state. We will therefore study first the simpler case of the neutron 
momentum distribution and then turn to the more difficult proton momentum distribution. 
In the following equations will mean the following J2i±PNi- And whenever we omit the 
index i it means that we implicitly sum over all three particles. 

1. Neutron in s He 

In Eq. (|3lD, the sum over is constrained by the energy conserving 5-function, and 
for the neutron spectrum in 3 He this gives a scattering state for the final two protons with 
the neutron off-shell. As a result the neutron does not satisfy the on-mass-shell relation 
E 2 = p 2 + m 2 . Since we are using a non-relativistic wave function for 3 He we will use a 
non-relativistic approximation for the relation between the energy and the momentum. We 
then define the binding energy of the nucleus, E, by the relation M = 3m + E, where m is 
the mass of a nucleon. Since we are working with a non-relativistic wave function, we make 
use of the approximation p° rs m + p 2 / (2m). As a result, the energy of the struck nucleon 
is p° a = m + E — p 2 p/(2m) — p 2 /(2m). One then finds p° a = m + E — p 2 t /(2jj) — q 2 a /(2v), 
where v is the reduced of the mass of the interacting pair and fi is their total mass 2 . If 



we compare this result with the expression given in Eq. (|3lD , then the recoil energy Tr is 
p 2 a / (2/i), while the separation energy, e, is E — q 2 a / (2z/). So the unpolarised spectral function 
for the neutron in 3 He is given by 

Snip) = \ E / °A (p, q) \Pn\ <TA (P, ?)) (35) 

® A 



2 Note that here, in the case of two identical particles we have v = m/2 and fi = 2m. 
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m + hi . 

2/i 2v) ) 

We stress that the two forms of Eq. (|3Tf ) are equivalent and should give the same results. 
In order to demonstrate this we computed the light cone momentum distribution, using 
Eq. © 

f n (y) = J d*k U + !pj8(y- S n (k) } (36) 



with the two forms of Eq. (pT|) . For the second form of this equation, the final state \<p, <r&) 
was taken to be a plane wave plus a pair of proton interacting in the 1 S'o channel. This is 
by far the most important channel for the final state interaction. We found that the light 
cone momentum distributions computed with the two forms of Eq. flST] ) were identical, for 
all purpose. 

For the polarised case there are two useful spectral functions 



#(P) 




(37) 



,6U- U + E-t-- h]] . (38) 



These spectral functions are, respectively, for a neutron with spin parallel or antiparallel to 
the spin of the nucleus. The "+" designates a positive projection of the spin of either the 
neutron or the nucleus on the z-axis, and the "— " a negative projection. In the same way 
as we obtain f n (y) we can calculate the quantities, f^(y) an d fn(y), just by inserting the 
correct spectral functions. Then one can form the useful quantity Af n (y) = f^(y) — fn{y), 
which is the equivalent of f n (y) for polarised structure functions. 



2. Proton in 3 He 



In the case of the proton we have two possibilities for the final state, so we also have two 
spectral functions. The first state is a scattering state similar to the final state encountered 
in the neutron case, with which it shares the formula for p°. The second possible final state 
is made of a scattered proton and a deuteron. We can find the form of the proton energy 
in the same way we did for the scattering state, only it is now much more simple as we 
have only two particles in the final state and not three. With the same non relativistic 
approximation as before, one easily finds that in this case: p° a = M — M d —p^/ (2M d ), where 
Md is the deuteron mass. Defining the binding energy of the deuteron, Ed, in same way we 
did for the tri-nucleon we have Md = 2m + Ed and finally, p° a = m + E — Ed — p 2 a / (2M d ). 
So we will have two spectral functions, Sp(p) (scattering state) and Sp(p) (deuteron state). 
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s p(p) = \ e/ & iflpi vI/± & «)) 




As in Eq. ( |3~7|) and Eq. (|38"D the "+" and "— " indicate the nuclear spin projection on the 
z-axis. 

In term of these spectral functions we can write the light cone momentum distribution 
of the proton 

f P (y) = lJ d4k ( 1 + ]p) 6 (y- ^pj ^ {k) + • (41) 

In the preceding equation we introduced a factor one-half because there are two protons in 
a 3 He nucleus. Without this coefficient f p would be normalised to 2 instead of 1. In the 
same way we did for the neutron we can extract polarised spectral functions, S^, for the 
proton by using a polarised density in combination with the right polarisation of the 
wave function. One can then get fp by applying Eq. fl4"T|), with the appropriate polarised 
spectral functions and in the end compute Af p (y) = f p {y) — f p (y)- 

C. Results 

Using the formalism presented above, we have computed light cone momentum distribu- 
tions for some of our three nucleon wave functions. For all those distributions we used only 
the first 42 three body channels. This because the computation of the polarised distributions 
involves some complicated matrix elements. However for all these wave functions the 42 first 
channels add up to more than 99% of the total, so one can safely assume that the contri- 
bution of the rest of the channels is negligible. For the unpolarised distribution the matrix 
elements are quite simple, so one can easily check, in this case, that the contribution from 
higher channels is indeed small. We compared the light cone momentum distribution for a 
proton and a neutron in 3 He for respectively 42 and 130 channels and found that for all pur- 
pose they were indistinguishable. For the PEST potential we also compared wave functions 
including 5 and 18 three-body channels and found that they were also indistinguishable. 
In Figs.[L] and |^ we show the proton and neutron light cone momentum distributions for 
our potentials (PEST, RSC and YAM7). The light cone momentum distributions given by 
the RSC and PEST potentials are almost indistinguishable and they cannot be separated 
on these figures. The YAM7 potential, however, shows some differences, probably because 
this potential does not include short range repulsion. It is also important to note that to 
have consistent results one needs to use a deuteron wave function computed with the same 
potential as the three nucleon system. 
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In Figs.|] and |] we show the proton and neutron polarised light cone momentum dis- 
tributions for the same potentials used in Figs.[l] and The polarised neutron light cone 
momentum distribution shows the same behavior and is similar in size to its unpolarised 
counterpart. However, for the proton the polarised momentum distribution is far smaller 
than its unpolarised counterpart. In this case all the potentials gives very similar results. 
We note that one can extract more information from the polarised momentum distributions. 
While in the unpolarised case the distributions are normalised to one, in the polarised case 
they are normalised to the polarisation of the given nucleon. From Ref. |2l| one can compute 
these polarisations analytically in terms of the S, S' and D waves probabilities (neglecting 
the small contribution of the P waves). One can compute those probabilities from the wave 
function and then compare them with the values extracted from the momentum distribu- 



tions. From Ref. p5| we have the following relations 



n + = I dyf+(y) = 1 - \ (P(S') + 2P(D)) , (42) 

n~ = I dyf~(y) = ± (P(S') + 2P(D)) , (43) 

P + = I dyf p + (y) = \-\ (P(D) - P(S')) , (44) 

P~ = fdyf-(y) = \ + \ (P(D) - P(S')) . (45) 

In Table |I| we compare the numerical values of these two expressions in 3 He, for our various 
potentials. The results in quite good agreement, with the small discrepancies arising from 
numerical errors in the computation of many nested integrals. (Note, for example, that the 



overall normalisation is correct to about 0.06%.) In Table |T| we make the same comparison 
but with wave functions in which we have adjusted the binding energies to the experimental 
values. 



IV. STRUCTURE FUNCTIONS 
A. Introduction 

In the incoherent impulse approximation, the structure function of a nucleus is the sum 
of the contributions from all its constituents. As we have already said in the previous section, 
the convolution formalism gives a way to link the in-medium structure functions to the free 
ones. This formalism, however has some limitations, especially at small Bjorken x, where 
other physics, like multiple scattering, becomes important. It is also only valid in the Bjorken 
limit, as the convolution formalism itself does not depend on Q 2 . In unpolarised scattering 
this formalism is a good tool to i nvestigate the EMC effect |26| , so we will use our previous 
results to study this effect in the three nucleon system. Another interesting result from the 
previous section is the fact that in 3 He, the proton polarisation (i.e. A p = p + — p~ « —2%) 
is very small and negative, while the neutron polarisation (i.e. A n = n + — n~ « 87%) 
is quite big. This is also clear from Figs.|3] and |j. This means that the neutron carries 
most of the spin of 3 He, so, at least for polarised scattering, this nucleus should be a good 
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approximation to a pure neutron target. The same argument is valid for the proton in 3 H. 
Since we already have a free proton target this may appear less interesting at first sight. On 
the other hand, it provides an ideal way to study the effect of the nuclear medium on the 
spin structure of a bound nucleon. 



B. Unpolarised structure function and EMC effect 

As we explained at the beginning of the previous section, in unpolarised deep inelastic 
scattering of a charged lepton on a nuclear target, all the target information is included in 
the two structure functions F\ and F 2 . In a simple quark model those functions have the 
following form [13,15] 

F 1 (x,Q 2 ) = W e2 M x >Q 2 )> ( 46 ) 
A q 

F 2 (x, Q 2 ) = 2xFt(x, Q 2 )=xY J e 2 q q(x, Q 2 ). (47) 

1 

In these expressions q(x) is the distribution of quarks of flavour q and electric charge e q . 
The relation between Fx and F 2 implies that the partons have spin 1/2 and no transverse 
momentum in the infinite momentum frame. A more general relation between F\ and F 2 

g is 



F 2 {x) = 2xF 1 (x) \ +R , (48) 
1 + 2xm N jv 

where R is the ratio of the cross section for absorbing a longitudinal photon to that for a 
transverse photon. 

Given the relation between F\ and F 2 , most studies concentrate on the latter. The 
convolution formula between the free and in medium F 2 structure functions [p^Pl is 

M A /„ \ 



Fi\x } Q 2 ) = dyf N (y)F» (j-,Q^| . (49) 

Hence the F 2 structure function of a nucleus of mass number A and proton number Z is 
given by 



F 2 A (x,Q 2 ) = I J dy (zf p (y)F$ ^,Q 2 ) + (A - Z)f n (y)F? (^Q 2 )) ■ (50) 



In comparing the F 2 structure functions on various targets, the European Muon Collab- 
oration (Aubert et al. f26|) discovered what is now called the "EMC" effect. We define a 
theoretical EMC ratio as the ratio of the F 2 structure function of the nucleus to the sum of 
the free structure functions of the nucleons in this nucleus: 

R t = F 2 A /(ZFi + (A-Z)F?) (51) 

On the other hand, it is more common to compare the ratio of the F 2 structure function of 
the nucleus to that of deuterium: 
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R x = (F 2 A /A)/(F 2 D /2) (52) 

This should be close to R t if the deuteron is a quasi-free system of a proton and a neutron 
and if the nucleus studied is symmetric, or almost, in its content of neutrons and protons. 
3 He and 3 H are highly asymmetric nuclei, as their content in one type of nucleon is twice as 
much as the other. To take this into account, it is common to an isosymmetric correction 
so that the ratio studied is 



F 2 D (x } Q 2 ) 



with 



Ur O 2 ) = ^,Q 2 ) + i?(x,Q 2 ) 

1[X ^ } ZF%(x, Q 2 ) + (A - Z)F?(x, Q 2 ) ' [ ^ } 

This ratio is, strictly speaking, the ratio of the EMC ratios of the nucleus A and the 
deuteron. Following the same kind of procedure used in the previous section, one can com- 
pute the light cone momentum distribution of a nucleon in the deuteron. To be consistent, 
this ratio has to be computed with the same interaction for both the three nucleon sys- 
tem and the deuteron. To compute Ra we used several parametrisations for the quark 
distributions: 



The parametrisation "CTEQ5" from the CTEQ collaboration |27| . The collaboration 
gives several parametrisations, but we mainly used the one called "leading order", 
and it will be the one used when we talk about the CTEQ5 parametrisation, unless 
explicitly stated otherwise. 



The "GRV" parametrisation from Gluck, Reya and Vogt p8 



• The "DOLA" parametrisation from Donnachie and Landshoff 

These distributions are usually given for quarks in a proton and in order to compute neutron 
structure functions we used charge symmetry^ |Kj . In Figs. |5| and |6| one can see the ratio R% 
for 3 He and 3 H, with the CTEQ5 parametrisation at Q 2 = lOGeV 2 , for the three potentials 
studied. In Fig.|^ we show R3 in 3 He for the PEST potential alone but for all three quark 
distributions (again at Q 2 = lOGeV 2 ). We also studied the effect of adjusting the binding 
energy as described at the end of the first section but did not include it in Figs.^| and |6| 
because it would have confused the plot. This adjustment of the binding energy caused a 
slightly deeper EMC effect in both 3 He and 3 H and also a slightly steeper increase at high 
x. 



3 With the exception of the DOLA distribution which gives proton and deuteron distributions. In 
this case we took the neutron as the difference between the deuteron and the proton. 



17 



C. Polarised structure functions 



If one does experiments with both a polarised lepton beam and a polarised spin 1/2 
nuclear target, one needs two more structure functions, g\ and g 2 - One can perform various 
measurements of cross sections with several polarisations in order to extract those two 
structure functions. They are smaller than F\ and F 2 and g 2 , in particular, is often neglected. 
As we indicated in the introduction, the figures for the effective polarisation of the nucleons 
in the three nucleon system seem to indicate that the contribution to the nuclear spin 
structure functions from the doubly represented nucleons is severely reduced. Thus, this 
system should be a good approximation to a pure single nucleon target. At leading order, 
§\ has the following form [ |ll] . |3T| ,|32| 

gi(x,Q 2 ) = W4 A( l^Q 2 )- (55) 
A g 

In Eq. (p5[) , Aq are the polarised quark distributions. They involve the difference between 
the distributions of quarks with the same and opposite helicity from that of the nucleon. It 
is much harder to find a simple parton interpretation for g 2 [L4j| . 

The convolution formula relating the free spin structure function to that in-medium is 
the following: 



g?(x, Q 2 ) = I MJ -MN{y)g? (-, Q 2 ) • (56) 

We computed the g\ structure function of 3 He using the same three potentials as for F 2 . 
The results from those potentials are sufficiently close that we will only use the results from 
the PEST potential hereafter. To compute g\ we mainly used the NLO "standard scenario" 
of Ref. | 5cl|| . We also studied the impact of the off-shell correction from Ref. |M] on g\. (The 



off-shell correction was calculated using a local density approximation and the quark meson 



coupling model |35j to estimate the change of the parton distributions in a bound nucleon.) 
In Fig||we show the following three curves at Q 2 = lOGeV 2 : xg\{x) for the free neutron, as 
well as xgi{x) for 3 He with and without the off-shell correction. As one can see, the three of 
them are close. The main complication in the extraction of g\ for the free neutron from 3 He 
is that the free proton spin structure function is very big compared with that of the neutron. 
So, while its contribution in 3 He is severely reduced by the low effective polarisation, it is 
still not negligible. One way to estimate the size of the contribution of the proton is to 
compare gi( 3 He) with a formula often used in the experimental analysis (see Ref. [15) for a 
derivation): 

9l ( 3 He) « A n9l (n) + 2A p9l (p). (57) 

If the contribution of the proton to gi( 3 He) is negligible, Eq. (|57f) is equivalent to: gi( 3 He) ~ 
A n g\{n). To estimate the effect of the proton contribution in the extraction of gi(n), we 
plotted the following differences: 

g 1 ( 3 He)-2A p g 1 {p) 
A <? = t 9i(n) (58) 
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and 



A g = — \ 9i{n). (59) 

In Figs.|] and |H] we plot both A g and A^. The second plot includes the off-shell effect of 
Ref. pi| . Note that the curves have been divided by / dxgi( 3 He) (~ — 1/16) so that one 



can judge the effect on the spin sum rule. Since one ultimately wants to extract <7i(n), we 
have also plotted that with the same normalisation, so as to have an idea of the size of the 
error in the differences^. It is clear from both plots that one gets more accurate results by 
including the proton contribution for mid-range x (0.2 < x < 0.6), the biggest error in this 
region occurring when the structure function crosses the x-axis. At higher x {x > 0.6) the 
effect of Fermi motion is significant and this will be even more important for 3 H, below. 
Nevertheless, the absolute value of the structure function is small and the corrections have 
little effect on the spin sum rule. For smaller x (x < 0.2) one clearly needs some other 
tools to accurately extract the free neutron structure function, even if the relative difference 
seems to be small. We find similar curves for other parton distributions, such as those from 



Ref. [36 



In the case of tritium one can plot a ratio, as gi(p) does not change sign. Therefore, to 
illustrate the effect of the neutron contribution in this case we plot: 

_ m ( 3 H) - 2A„ gl (n) 

R » a^IF) (60) 



and 



9 A P 9i(p) 



In Fig. |TT| we show both ratios (R g is the solid line and R g is the dashed line) without 
including the off-shell corrections as well as R g with the off-shell corrections (dot-dashed 
line). In this figure we can clearly see that on most of the interval the contribution of the 
neutron is negligible, some difference appearing for small x. This is expected simply because 
gi(n) is significantly smaller than gi(p) for most values of x. On the other hand, we can 
also see that medium effects seem to be quite important and that the off-shell correction 
makes an important difference. One can also see clearly the effect of Fermi motion at high 
x, while it would be invisible if one were to plot differences. It is clear from these results 
that from a measurement of gi( 3 H) one can expect to extract the size of the change in the 
spin structure function of the bound proton and one might even hope to separate the origin 
of this effect. 



4 We do not plot the ratio of structure functions because in both the neutron and 3 He cases g\ 
can be zero, leading to singularities in the plots. 
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V. CONCLUSIONS 



We have computed the three-nucleon structure functions from various two body poten- 
tials. This involved calculating wave functions, light cone momentum distributions and 
finally the structure functions. We have presented our prediction for the EMC effect in both 
3 He and 3 H and shown that they were quite close for various two-body potentials and quark 
distributions. In addition, we saw that isospin breaking would have only a small effect on 
these findings. This result has been used elsewhere [57| in a proposal to measure the d/u 
ratio at large x at Jefferson Laboratory p8| , |39 . 



From our study of the spin structure function of 3 He, we showed that it possible to 
extract the structure function of a polarised neutron with reasonable accuracy. However, it 
is necessary to account for the contribution from the pair of protons which are not totally 
unpolarised. Turning to the polarised structure function of 3 H, we saw that while the exper- 
iment is extremely challenging it could also be very valuable. In particular, one can measure 
the size of the medium corrections and check experimentally the predicted modification of 
the spin dependent parton distributions of the bound nucleon. 
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APPENDIX A: THE KERNEL OF THE HOMOGENEOUS FADDEEV 

EQUATION 

For completeness, we present in this Appendix the explicit expression for the kernel of the 
homogeneous Faddeev equation when the interaction is represented by a separable potential. 
The details of the derivation are in Ref. [01. We have 



7JI - I n ria. r>J! 



dx ~ FT! f V T e a N a ;e g N s {Pa,Pp;X) , (Ai) 



where 



1 f dx 9%(q«)9%k 

2 l X i (j'n ■ Pfi ■ I' ■ 

r £N a -ApNp (p«> p& x ) = ( ? J ( ~r ) Bn « n ? ^ Pc ^ 



.9a j \qp. 



C 




xVV A c ^ b ^ , (A2) 

a=0 b=0 \PP , 



with Pc(x) the Legendre polynomial of order C, and 

X=Pa-P(3 qa = -Pp ~ ^Pa qp=P a + 2^/3 • ( A3 ) 
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The coefficients A^' ^ t N which results from the recoupling of the spin and orbital angular 
momentum is given by 



.4 



C,a,b 



l) R £ a £pL a Lp S a Sp j a jp s a §p C' 2 p a a pp 



2 „a J> 



(24 + 1)1(2^ + 1)! 







!) 
















A 





\J (2a)! (26)! (24 - 2a)! (24 - 26)! 
x E (/AA') 

/AA' 

ja S a Sp jp 

Ja f Jf3 h 

S a 4 4 S /3 

4 4 / 
a 4 - b A 
6 A' 



x 



X 



X 



A' C Lp \ ( A £ L a 



a 4 - 6 A 




, - a 6 A' 




where the 12 — j symbol is that defined by Ord-Smith | i0|| , the phase R is defined as 

R=-J + L a + Lp + S a + Sp+j a +jp-j a + sp + £ a + C, 
and finally p a and pp are 

1 m, 



(A4) 



Pa 



P/3 



1 

2' 



rag + m 7 2 ' ' K ra a + m 7 
The isospin recoupling coefficient B^ a ^ f} is given in terms of 6 — j symbol by the relation 



B 



NaNp 



2/3 ^7 la 
l a I lp 



(A5) 
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TABLES 



TABLE I. binding energy for a given potential and components of the wave function. 



Potential 


number of 


binding energy 


P(S) 


P(S') 


P(D) 




channels 


(MeV) 


% 


% 


% 


RSC 


5 


-7.15 


88.37% 


1.88% 


8.89% 


YAM4 


5 


-9.12 


93.08% 


1.58% 


4.97% 


YAM7 


5 


-8.05 


89.1% 


1.59% 


8.71% 


PEST 


5 


-7.27 


89.3% 


1.88% 


8.11% 


PEST 


10 


-7.10 


89.72% 


1.71% 


7.85% 


PEST 


18 


-7.32 


89.56% 


1.66% 


8.07% 



TABLE II. Effective polarisation of the nucleons in 3 He for various potentials. 













If(y) 






n + 


n 


p+ 


p- 


n + 


n 


p+ 


p- 


PEST 


93.97% 


6.03% 


48.96% 


51.04% 


93.62% 


6.32% 


48.98% 


50.96% 


RSC 


93.45% 


6.55% 


48.83% 


51.17% 


92.92% 


6.79% 


48.76% 


50.95% 


YAM7 


93.66% 


6.34% 


48.81% 


51.19% 


93.25% 


6.35% 


48.69% 


50.92% 



TABLE III. Effective polarisation of the nucleons in 3 He and 3 H, with two-body interaction 
adjusted to produce the experimental binding energies. 













If(y) 






n + 


n~ p + 


p~ 


n+ 


n p + 


p~ 


3 He 


93.97% 


6.03% 48.91% 


51.09% 


93.73% 


6.24% 48.94% 


51.02% 


3 H 


93.45% 


6.55% 48.85% 


51.15% 


93.86% 


6.13% 48.89% 


51.10% 
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FIGURES 
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FIG. 1. Neutron light cone momentum distribution in 3 He for various potentials. 
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FIG. 2. Proton light cone momentum distribution in 3 He for various potentials. 
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FIG. 3. Neutron polarised light cone momentum distribution in 3 He for various potentials. 
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FIG. 4. Proton polarised light cone momentum distribution in 3 He for various potentials. 
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FIG. 5. The ratio i?3, given in Eq.(|53|), for 3 He, at Q 2 = lOGeV 2 , calculated for various 
potentials using the CTEQ5 quark distributions. 
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FIG. 7. The ratio R 3 , given in Eq.flU, for 3 He, at Q 2 = lOGeV 2 , calculated for the PEST 
potential, using various quark distributions for the nucleons. 
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FIG. 8. Comparison of several calculations of xg\{x) for 3 He, at Q 2 = WGeV 2 . 



28 





FIG. 10. A g , A' g and gi(n), including off-shell corrections, at Q 2 = WGeV 2 . Note that all 
three curves have been divided by / dxg\( 3 He). 
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